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Abstract 



o 

Q 

Q> I We compute the photon production rate of a quark-gluon plasma (QGP) 



^ at finite quark chemical potential /i using the Braaten-Pisarski method, thus 

' continuing the work of Kapusta, Lichard, and Seibert who did the calculation 

o ' 

■ for ;U = 0. 



The thermal production of hard photons in ultrarelativistic heavy ion colhsions has been 
proposed as a possible signature for the formation of a quark-gluon plasma (QGP) Since 
the mean free path of photons in the fireball is much larger than its dimensions 0, photons 
provide a direct probe of the fireball. So far the production rate of hard photons has been 
considered at finite temperature but vanishing quark chemical potential mainly. In order 
;h ' to include medium effects consistently the Braaten-Pisarski resummation technique has 
been applied to this problem 

Assuming the formation of a QGP already at AGS and SPS energies, however, a finite 
quark chemical potential has to be considered 0, and even at RHIC energies the quark 
chemical potential /i may not be negligible as indicated by RQMD simulation \ ~2T) 
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0. At given energy density estimates based on lowest order perturbation theory indicate 
a strong suppression of the photon production at non-vanishing fi as compared to the case 
/i = 0|. The aim of the present work is the improvement of these estimates by applying 
the Braaten-Pisarski method generalized to finite chemical potential 0. 

We follow the calculation of Kapusta, Lichard, and Seibert [Q, and describe the changes 
that have to be made for a nonzero quark chemical potential. For brevity, we do not 
demonstrate the intermediate steps of the computation where they would be analogous to 
ref. i. 

To leading order the hard photon production rate is derived by using the Braaten- 



Yuan prescription [|T0| resulting in a decomposition into a soft part, which is treated using 
the resummed propagators of Braaten and Pisarski, and a hard part containing only bare 
propagators. For this purpose a parameter kc is introduced, separating the soft from the 
hard momenta of the intermediate quark. Demanding gT <^kc'^T and assuming the weak 
coupling limit, the final result is independent of the separation scale kc- 

The soft part can be obtained from the imaginary part of the self energy of a photon 
propagating through a QGP. The contributing diagrams to the photon self energy are shown 



in fig. 1. There quark lines with blobs represent effective quark propagators [Til], however 
for the /i 7^ 0-case with a modified quark mass [^|12| of 

,2 / ,,2\ 



Our result for the soft contribution to the production rate of a hard photon with energy E 
and momentum p is 

AR'"'' 2 „„. 1 



2E—— = -- — -^g'^'Ini U 



This differs from the result obtained in ref. only in that (a) the quark mass nig depends 
now on the chemical potential and (b) the factor in front of the logarithm is modified 



since it originates from a factor m^. (In deriving eq. (|^) we assumed a hyperbolic cut-off 
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> k — CO > in accordance to the hard part (y) following ref. where u is the energy 
and k the momentum of the spacelike intermediate soft quark. Thus uj and \k\ need not to 
be soft individually. However, closer investigations show that for the expression under the 
loop integral in H^^*, approximations are legitimate that assume u as well as \k\ to be soft, 
since it is this region where the main contribution to the loop integral comes from.) 

The easiest way to calculate the hard part is starting from the scattering matrix elements 
shown in fig. 2. They are related to the photon self energy of fig. 1 by cutting the latter 
13| . From QCD Feynman rules, we find the squared matrix elements to be 



= (3) 

for the annihilation process, and 

= _^vr2aa,^^±^ (4) 
9 st 

for each of the two Compton processes for quarks and antiquarks. The symbol S indicates 
that these matrix elements are already summed over spins, colors, and two flavors {u and 
d). The letters s, t, and u denote the Mandelstam variables. In each diagram, all legs but 
the outgoing photon belong to thermalized particles. The photon is not thermalized since 
it has a large mean free path. For each of the three possible processes we may compute the 
hard part of the photon production rate as |Q 

(5) 

where ^11,2,3 are Bose or Fermi distribution functions, respectively; the plus sign is for the 
annihilation process and the minus for the two Compton processes. Still following ref. [Q, 
we may rewrite this equation as 
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with the polynomial P{Ei, E2) = -itEi + {s + t)E2f + 2Es{{s + t)E2-tEi)-s'^E'^ + sH + st'^ 
and the step function O. 

For the /i = 0-case, it is a reasonable approximation to use Boltzmann distribution 
functions for rii and n2 instead of the full Fermi/Bose functions. In this way, all integrations 
may be performed analytically, and the result 

dR 5 aa.T^e"^/^ / 2 , fAET\ 1 f AET\ \ 

V ' " V ' 

annihilation Compton 

fits perfectly together with the soft part: the sum of both is independent of kc- 

Actually, the Boltzmann approximation (for ni and 712) is better than it should be, 
considered that the photon contribution from the Compton processes is significantly un- 
derestimated (up to 30%) and the annihilation contribution is overestimated. Surprisingly, 
both errors cancel each other, up to an error of about 10% in the final result for those values 
of T, E, and g that are interesting in practice. 

The total photon production rate for = 0, computed in the Boltzmann approximation 
(for ni, n2), is 

which also follows from the photon damping rate by the principle of detailed balance [Q. 

Dumitru et al. use the Boltzmann approximation also in their computation of the 
photon rate at /i 7^ 0, however they obtain a hard part which does not match onto the ana- 
lytically known soft part, (||). The cited result also contains a term ~ Ei ((4i?/i — kD/AET) 
which runs over a branch point at finite /i. 

Our approach to the problem is the following: since we cannot evaluate the integrals in 
(^) containing the exact distribution functions analytically, we employ numerical methods. 
We rewrite (|^) in a form suitable for Gauss quadrature (i?+ := Ei + E2): 
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Laguerre Legendre 
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zp g-iE+-E-f,3)/T 



(9) 
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Laguerre 

dE2 



g-Mi/T _|_ g-(E+-E2)/T g-/i2/T _p g-B2/r 



Chebyshev 

In this compact notation, the upper signs and /ii — —^2 — /J', /J'S — apply to the annihi- 
lation process; the Compton processes require the lower signs and fii = 0, /X2 = = ±A*, 
where the two results for +11 and —/i have to be added in order to take antiquarks as well 
as quarks into consideration. The polynomial Pi(£'+, E2) is just P(-E'+ — -E'2, E2)/s^] thus it 
has the leading coefficient —1 and gives under the square root a perfect weight for a Gauss- 
Chebyshev quadrature in E2. The £^+-integral, as well as the s-integral, is done numerically 
by a Gauss-Laguerre quadrature; the necessary exponential weight function arises naturally 
in the expression. The i-integral, having no appropriate weight function, is performed via 
Gauss- Legendre quadrature. Finally, the s-integral is done by a Gauss-Laguerre quadrature, 
as suggested again by the naturally arising weight function. 

Unfortunately, each of the four integrands is highly peaked at the ends of the integration 
interval, due to singularities in or near the domain of integration. This means the inte- 
gral value is dominated by contributions coming from comparably small regions, and the 
quadrature problem is ill-conditioned. In order to cure the problem, we subtract from the 
integrand in each step of the fourfold integration those contributions that stem from poles 
inside or slightly outside the integration region. Those parts are integrated out separately 
in an analytical way; the remaining numerical integrals are much easier to compute, since 
the integrands are no longer varying strongly. This way, we need only about 20 points for 
each quadrature, and still obtain numerical results for the hard part that are precise up to 
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an maximal error of 0.2%. 

In fig. 3 kc is varied over almost six orders of magnitude, leaving the total photon rate 
fixed within numerical error. (In order to test the cancellation of kc, expected to occur in 
the weak coupling limit, we adopt a value of g = 0.01.) Even in regions where gT <^ kc '^T 
is by no means true, the sum of soft and hard part is perfectly constant, while the soft part 
alone varies over a wide range and even becomes negative. This is a nice verification of the 
Braaten-Yuan method for 7^ P,|lOl . From (|^) and the independence on kc we know that 



the final result has to assume the form 



where the dimensionless quantity G follows from the hard part. As a result of the numerical 
calculation of the hard part, discussed above, it turns out that G depends only on fi/T. 
Within a 3% error in the rate, G depends very weakly on E/T. However, G is nicely 
independent of T, as a dimensionless quantity should be. We found that G can be fitted to 
a good approximation by the phenomenological formula G = ln(l + yU^/vr^T^), thus leading 
to the compact expression 

This pocket formula reproduces the exact, numerically calculated result within an error 
of 3% for |/i/T| < 1. A larger fi/T requires a one-parameter fit for G. We found that 
G = ln(l + 0.139;U^/T^), inserted in eq. ([T0|), leads to a phenomenological formula for the 
rate that is precise on the 3%-level even for much larger |/i/T|. We emphasize that the error 
of 3% is not due to the numerical evaluation of the hard part but to the assumption that G 
is solely dependent on fi/T. 

We are now able to extrapolate this formula to realistic values of g and discuss the 
photon spectrum of the QGP. For photon energies E above about 3T the logarithm in ( pA]) 
is positive, indicating the validity of the extrapolation to realistic values of the coupling 



constant |Tj] . Fig. 4 shows how the spectra are dominated by the exponential decrease with 



temperature. We draw the conclusion that a measurement of the photon spectrum may serve 
only or mainly to find out the temperature of the plasma; in order to measure the chemical 
potential by a photon experiment, one would have to measure the overall coefficient of the 
exponential, which is in an experimentally obtained spectrum infiuenced by many other 
factors not yet under control, like the size and duration of the plasma phase. However, if 
one has information on the energy density e of the plasma, T and fi are related to each other 
by an equation of state like the one quoted in ref. 



V 30 3 J V TT y ^ 27r2 V TT 

If T is made dependent on fi in this fashion [§, the resulting photon spectra are strongly 
dependent on fi. At RHIC, one expects a maximum energy density of about e = 5 GeV/ fm?. 
In fig. 5, we show the corresponding photon spectra. The temperatures of the five curves 
decrease with rising /i from Gd^ to about 0.22 Gel^, whereas n/T varies from to about 
1.8. The photon suppression at finite chemical potential and fixed energy density observed 
in ref. ^ is therefore an indirect phenomenon, caused by the reduction of the temperature. 

In conclusion, at /i = 0, our result ( pA]) conincides with the one in ref. |Q. The Boltzmann 
approximation used in there seems to work fine. 

At /i 7^ 0, Bose and Fermi distributions have to be used in order to match the soft 
contribution onto the hard one according to the Braaten-Yuan prescription providing a 
consistent result to leading order. Hence the hard part of the photon production rate 
can only be determined numerically. We were able to demonstrate that the final result is 
independent of the arbitrary separation parameter kc and obtained the simple formula (pTD. 
Our graphs of the photon spectrum have a strong similarity to the one of Dumitru et al. , 
mainly because they all are dominated by the single factor e~^l^ . 

We did not consider pre-equilibrium effects so far. The equilibrium distribution func- 
tions used in our calculation are not quite applicable, since the plasma is probably never 
completely equilibrated. For a phenomenological treatment of a pre-equilibrium plasma, 
fugacity factors may be employed JlSl . This will be dealt with in a future publication. 
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